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The surface of a 3+ld topological insulator hosts an odd number of gapless Dirac 
fermions when charge conjugation and time-reversal symmetries are preserved. Viewed as 
a purely 2+ld system, this surface theory would necessarily explicitly break parity and 
time-reversal when coupled to a fluctuating gauge field. Here we explain why such a state 
can exist on the boundary of a 3+ld system without breaking these symmetries, even if 
the number of boundary components is odd. This is accomplished from two complemen- 
tary perspectives: topological quantization conditions and regularization. We first discuss 
the conditions under which (continuous) large gauge transformations may exist when the 
theory lives on a boundary of a higher-dimensional spacetime. Next, we show how the 
higher-dimensional bulk theory is essential in providing a parity-invariant regularization 
of the theory living on the lower-dimensional boundary or defect. 
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1. Introduction 

It is well known that a 2+ld theory consisting of an odd number of gapless Dirac 
fermions interacting with a fluctuating gauge field must break parity (P) and time-reversal 
(T) symmetries [|l],l|,EI • It is also well-established that there exist topologically non-trivial 
band structures for fermions on a lattice in three spatial dimensions (i.e., 3+ld topological 
insulators) whose surfaces harbor an odd number of gapless Dirac fermions [dlUJ. This 
raises the question: Must topological insulators break P and T on their boundaries when 
coupled to a fluctuating gauge field? The purpose of this paper is to explain why this does 
not occur. 

Strong topological insulators in 3+ld (3DTI) are distinguished by the presence of 
gapless surface states that are stable to all P and T invariant perturbations (with respect 
to the boundary theory) that conserve electric charge @. Provided that the chemical 
potential is fine-tuned to the Dirac point, these gapless surface modes can be described at 
low energies by an odd number of 2-component Dirac fermions, which are charged under 
the electromagnetic {7(1)em § au S e field. In a purely 2+ld theory of this type (i.e., 
QED3) , gauge invariance (under both large and small gauge transformations) is preserved 
if and only if 

^+fcGZ, (1.1) 

where Nf is the number of flavors of 2-component Dirac fermions and k is the level of the 
Chern-Simons (CS) term for the gauge field, i When Nf is odd, a non-zero half-integral 
level CS term must supplement the effective action. 

This is known as the parity anomaly. It is the gauge-invariant regularization of the 
theory that results in the addition of the half-integral level CS term to the effective action 
when Nf is odd. This CS term explicitly breaks P and T. Physically, the anomaly means 
that parity and time-reversal invariance are explicitly broken when a theory with an odd 
number of 2+ld Dirac fermions is coupled to a fluctuating gauge field. 

1 To be specific, the U (1)em gauge field is not confined to the 2+ld surface; it is free to explore 
the 3+ld bulk. Integrating over this extra direction results in a tree- level propagator for the gauge 
field that is softer in the infrared (IR) in that it diverges linearly as opposed to quadratically at 
small momentum. 

2 This statement is slightly imprecise. Invariance of the theory under large gauge transforma- 
tions requires a choice of auxiliary 4-manifold to ensure the CS term is well defined in general 
[Q^] . While this more precise definition is suggestive of our eventual conclusion, we will not make 
use of it further. 
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Anomaly considerations establish a relationship between certain topologically ordered 
phases of matter and the gapless modes living on their boundaries |T0| , Pi|T2| , p!3[1 . The quan- 
tum Hall effect provides, perhaps, the most famous example of how anomaly considerations 
can be used to better understand these gapless boundary modes in a model independent 



way |TT|JT^]. In this example, the £7(1) em charge conservation symmetry is gauged by 
the electromagnetic field, whose (possibly fractionally) quantized low energy Hall response 
implies an effective bulk description by a CS theory. Equivalently, the non-dissipative Hall 
current requires charge-carrying chiral edge modes. The potentially anomalous £7(1)em 
gauge symmetry provides the link between these two descriptions. In the presence of a 
boundary, neither the bulk CS theory nor the boundary chiral theory is individually gauge 
invariant; however, their anomalous variations cancel one another so that the underlying 
£^(1)em symmetry is maintained in the system as a whole. Thus in the presence of a 
boundary, the bulk theory cannot exist without the gapless boundary modes, and vice 
versa. 



This effect is known as anomaly inflow ||14|1 . A classical anomaly of the bulk effec- 
tive action is cancelled by a local quantum mechanical anomaly of the boundary theory. 
Intuitively, the cancellation occurs because charge flows out from the bulk and along the 
boundary of the system at a rate determined by the applied external field. Because this 
relationship readily extends to interacting systems, it is a powerful demonstration of the 



robustness of the gapless boundary modes [13 



Anomaly inflow does not, however, define the relationship between bulk and bound- 
ary theories for 3+ld topological insulators. This is because local violation of charge or 
momentum conservation can only occur when the spatial dimension is odd [|IT^|l6l , |T7| ,j3f , 
as is the case for the 1+ld boundary of a 2+ Id quantum Hall system. Instead, possible 
anomalies relevant to theories in even spatial dimension necessarily involve so-called large 
gauge transformations. This is the case for the parity anomaly constraint ( |1 . 1|) , which is 
potentially relevant to the surface of a 3+ld topological insulator, i 

A (topological) band insulator is necessarily a system that can be realized on the lattice 
so it is worth reviewing the conventional wisdom regarding anomalies in lattice systems. 
Because the lattice itself provides a gauge-invariant regularization, any system that can 
be realized on the lattice cannot be anomalous. That this is true follows immediately 



from fermion doubling [20] . For example, a purely 2+ld lattice system of fermions with 



See [18,19] for an earlier discussion on a closely related system. 
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relativistic dispersion necessarily contains an even number of (low energy) Dirac fermions 
in the absence of P and T breaking. The parity anomaly constraint ( |1 . 1|) is satisfied by 
the low energy effective theory for all deformations preserving the U(l) symmetry - even 
those that break P or T ||21|| . The striking feature of topological insulators is that the 
gapless surface modes do not exhibit fermion doubling because they live on the boundary 
of a higher-dimensional system. 

Thus, we must ask whether ( |1 . 1| ) is obeyed on each boundary or, if not, how this is 
consistent with an overall gauge-invariant theory. Because of the nature of the parity 
anomaly, a topological insulator with two (separate) boundaries automatically satisfies 
( |1.1|) . However, there is no fundamental reason why a topological insulator must have 
an even number of boundaries; for example, a solid sphere or torus of 3+1 d topological 
insulator would have a single boundary, potentially violating (|1.1|). 

In this work, we explain why ( |1 . 1[ ) need not be satisfied by the low energy theory 
describing surface Dirac fermions interacting with a bulk gauge field. Specifically, we shall 
explain how the higher-dimensional bulk theory from which the surface modes descend 
eliminates the potential anomaly of the surface. We emphasize that this conclusion is 
quite different from what occurs in the quantum Hall case [O] and other examples studied 



by [13], where the potential anomaly of the surface is cancelled by a comparable anomaly 
in the bulk. 

In fact, we find that ^y- + k is half-integral at each boundary surface of a 3+ld 
topological insulator only in the limit that the bulk gap mo is infinitely large compared 
to any T-breaking perturbations on the boundary. The corrections appearing at finite mo 
imply that the CS level k need not be quantized at half-integral (or indeed at any rational) 
values, even in a non-interacting system. 

The conclusion that the surface states do not exhibit the parity anomaly -in the 
sense that they do not obey ( |1 . 1| ) - is of clear importance for the low energy properties 
of a topological insulator: a non-zero bare CS term breaks P and T, and would therefore 
drastically affect the low energy physics. While these surface properties are well established 
(or at least well believed) theoretically for the case of topological band insulators which 
need not be coupled to a fluctuating gauge field in order to be defined, our analysis is 
equally applicable to the case of more exotic fractional topological insulators, such as 
those described in p^ , p3| , p3] , |25|1 , in which the presence of a fluctuating "internal" gauge 
field is inevitable. Thus, our result is important in clarifying why P and T invariant gapless 
boundary modes exist in these systems - as well as in understanding their topological order. 
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We would like to point out that in the context of topological insulators, a different 
definition of the parity anomaly is sometimes used. It is sometimes said that the surfaces 
modes exhibit the parity anomaly because the parity-violating current equation, 

W = ^pPA", (1-2) 

is satisfied on any boundary perturbed by a P and T odd interaction. For a topological 
insulator boundary, ( |1.2| ) does not mean that there exists a zero magnetic field Hall con- 
ductance (in contrast to a purely 2+ld system); rather, a Hall effect occurs only after a 
time-reversal breaking perturbation has been applied to the surface, as we shall explain. 
( (|1 . 2| ) follows directly from the effective action calculated in Sec. 4.) In the present work, 
when we say that the surface modes of a 3+ld topological insulator do not exhibit the 
parity anomaly, we mean that they do not satisfy the constraint (|1.1|) requiring integral 
Nf/2 + k (on each boundary component). It is important to note this difference in termi- 
nology in order to avoid possible confusion. @ 

The remainder of the paper is organized as follows. We begin in Sec. 2 with a brief 
introduction to the specific model we wish to study. We then turn to the explanation for 
why non-integral + k is consistent with a gauge-invariant theory. To do so, we make 
use of two complementary perspectives: topological quantization conditions in Sec. 3 and 
perturbative regularization in Sec. 4. We summarize and conclude in Sec. 5. 

The paper contains five (count them!) appendices summarizing issues that are related, 
but not essential to the above line of argument, although they may be of some interest. 
Appendix A contains a review of domain wall fermions and their relation to continuum 
models of topological insulators. In Appendix B, we discuss anomaly inflow intuition for 
line and domain wall defects in three spatial dimensions. In Appendix C, we recall how flux 
insertion arguments can be used to define a strong topological insulator in the presence of 
disorder or other interactions. In Appendix D, we repeat the perturbative analysis of Sec. 4 
in the technically simpler, but conceptually equivalent case of 1+ld allowing direct contact 
with the work of Goldstone and Wilczek |2(J . In Appendix E, we elaborate in detail upon 



the leading divergence structure arising from the interaction between bulk and boundary 
modes in our toy model of a topological insulator, thereby confirming the conclusions of 
Sec. 4. 



We thank E. Fradkin for discussions on this point. 
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2. Domain Wall Fermions, the #-term, and Topological Insulators 

There are four Z2 topological invariants that characterize the bulk band structure 
of (non-interacting) fermionic insulators in 3+ld [[|,[5|,[| . These distinguish between three 
classes of time-reversal invariant, charge- conserving band insulators: a so-called trivial 
insulator (with no protected low-energy surface modes), a "weak" topological insulator 
(which has gapless surface states that can be gapped without breaking time-reversal sym- 
metry, but are nonetheless robust to disorder [p7p, and a "strong" topological insulator 



(STI), whose gapless surface states cannot be eliminated by any time-reversal invariant 
perturbation, i In this paper, we focus on the STI. 

At low energies, a STI in 3+ld can be described by a continuum theory of a single, 
massive 4-component Dirac fermion [ FZoySOfl . In this continuum formulation, the STI is 
distinguished from its trivial counterpart by the sign of the Dirac fermion mass m. The 
existence of two distinct insulators distinguished by the sign of m is the continuum version 
of the notion of topological band structure: these two insulators cannot be adiabatically 
connected without either closing the bulk gap or choosing a connecting path in parameter 
space that breaks P and T (e.g., by interpolating over complex fermion masses). 

To exhibit the difference, we consider the action, 



S = J (Tx^i^id^-ieA^ + mix^J, (2.1) 

where ip is a 4-component spinor that describes the massive bulk fermion with spatially 
dependent mass m(x) coupled to the U(l) gauge field A^ which may represent the elec- 
tromagnetic field. The spatially- varying mass allows us to study the interface between a 
topologically non-trivial and a topologically trivial insulator, where the parity anomaly 
constraint (|1.1|) could potentially be applied. (We have not included kinetic terms for 
the gauge field, since their specific form does not affect our results.) In (|2.1|), tp = "i/^7 , 
{7^7^} = for fx, v = 0, 1,2,3 where rf 11 = diag(l, —1, —1, —1). Further, we define 
the matrix 75 = Z7°7 1 7 2 7 3 which anti-commutes with all 7 M . 



5 It is well established that this is the case for non-interacting systems; however, the arguments 
of [28 1, applied in 3D as explained in [24], strongly suggest that this remains true in the presence 



of interactions. 
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Let us begin by reviewing the claim that the sign of the fermion mass distinguishes 
between two distinct bulk phases of matter. If m(x) is constant, the two states are distin- 
guished by the presence or absence of a topological #-term in their low-energy effective ac- 
tion [pO| , |3T|j32]1 . We can see this by considering the effect of chiral rotations ip — » exp^crys)?/? 
on the effective action ( |2.1|) . These rotate the fermion mass according to 



mipip — >■ mipe' l2ai5 ip = m cos (2a) -i/^ + im sin (2a) -075 "0 (2-2) 

rendering it complex unless a is an integer multiple of tt/2. Importantly, chiral rotations 
also contribute an anomalous term to the action from the path integral measure: 



e 
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S e = (2a) / d x—^e^p^F^Fp^. (2.3) 



327T 2 ^ upa ^ 9<J ' 

By the chiral rotation 2a = n - which is nothing more than a change of path integration 
variables - we may change the sign of the fermion mass at the expense of creating a 
topological #-term ( |2.3| ) with coefficient 9 = n. Thus, the effective continuum actions 
for the topologically non-trivial and trivial insulators with constant masses everywhere in 
space, but with opposite sign, differ precisely by the topological term ( |2~^ ) with 9 = tv. 

Now consider the scenario in which there is a single domain wall separating a region 
of STI (X3 > 0) from a region of the vacuum or trivial insulator (x^ < 0) . As we review in 
Appendix A, if we take 

lim m(xs) = ±mo, m > (2-4) 

X3 — >±oo 

with m(xs) passing through zero exactly once, at X3 = 0, there is a single massless 2+ld 
Dirac fermion localized near X3 = 0, where m(x) changes sign ||33|| . 

Based on the result ( |2.3| ), we might expect a #-term with a spatially varying coefficient 
taking the value, say, 9 = for X3 large and negative, and 9 = n for X3 large and positive. 
A bulk #-term integrates by parts to a boundary CS term at level k = 9/2n. This suggests 
that the low energy action might obey ( |l.l| )with the half-integer CS level, obtained from 
a bulk #-term, compensating for the odd number of domain wall fermions. If true, this 
would mean that P and T are explicitly broken at the surface of a topological insulator 
in the absence of any such symmetry-violating perturbations. In particular, a magnetic 
perturbation that opens up a gap in the surface fermion spectrum would imply an integral 
Hall effect, along with the consequent change of the Kerr and Faraday angles for light 



passing through a single surface pOJ33,R5 
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In the remainder of the paper, we show that this is not the case: we should think of the 
domain wall as associated either with the presence of an odd number of gapless 2+ld Dirac 
fermions, or (if these are gapped) with a half-integral surface Hall conductivity which may 
be understood as arising from a spatially varying #-term. This conclusion is equally true 
for theories of fermions coupled to other, possibly non-Abelian, gauge fields, as is relevant 
for describing strongly interacting fractional topological insulators. In other words, the 
surface of the topological insulator fails to obey ( |1.1|) ; + k is half-integral (in the limit 
mo —> oo). 



3. Large Gauge Transformations, the Parity Anomaly, and Topological Insu- 
lators 



In order to understand the applicability of constraint ( |1 . 1|) to the surface modes of 
a topological insulator, we first review its derivation from the perspective of topological 
quantization conditions in 2+ld. The generalization of this logic to topological insulators 
is then immediate. 



3.1. The Fermion Determinant 

Recall that there are two types of possible gauge anomalies of an effective action: local 
or global. (For general discussions, see PB|J^,P7|.) The distinction arises from the class of 
the particular gauge transformation under which the action fails to be invariant. If a local 
anomaly is present, the action fails to be invariant under any gauge transformation that 
is continuously deformable to the constant map. The current associated with a locally 
anomalous symmetry fails to be conserved. 

In contrast, the parity anomaly is an example of a global anomaly: an anomaly 
associated with so-called "large" gauge transformations. By a large gauge transformation, 
we mean a gauge transformation that is not continuously deformable to the constant map. 
For example, two maps with distinct winding number from the circle to itself cannot be 
continuously deformed into one another. 

Invariance of the effective action of QED3 under large gauge transformations requires 
that + k be integral. To see this, consider first the fermionic contribution to the gauge 
field effective action, 

[#][#] exp ( / tjj(iD { 3 A) - mo)V>) = (det(D { 3 A) - m )J ' , (3.1) 
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where = ^ a (d a —ieA a ) is the 2+ld Dirac operator for a given gauge field configuration 
A. The fermion determinant is essentially a product over the energies of the fermionic 
states. These energies are negative for states that lie below the chemical potential, and 
positive for states above it. (In this paper, the chemical potential is fine-tuned to zero so 
that charge-conjugation symmetry in maintained.) The square root instructs us to only 
include, say, the filled negative energy states in the product (p.l|). 

Therefore, an anomalous transformation of the fermion determinant (or anomaly, for 
short) occurs when an odd number of fermions are "pumped" from immediately below the 
Fermi surface to states immediately above it under a large gauge transformation. (Gauge 
transformations deformable to the identity can have no such effect.) In such a situation, 
the fermion determinant changes sign. This renders the partition function, which is a sum 
over all such sectors, ill defined if there is no compensating bare CS term. 

As a concrete example for how this works, consider a Dirac fermion on a spatial torus 
S 1 x S 1 . We will also assume that all gauge fields tend to constant values and so are pure 
gauge as t — > ±00, effectively imposing periodic boundary conditions in time on all physical 
observables. Suppose the fermions are given anti-periodic boundary conditions along the 
two cycles of length L\ and L 2 . Then the allowed fermion momenta (in the absence of any 
external gauge field) are k\ = f^-(^i + 2 )i k 2 = f^( n 2 + \) with associated band energies 
E = ±yfci~-j-~fcf. Consider the Dirac cone sitting at the time-reversal invariant momentum 
point (ki, k 2 ) = (0, 0). The allowed momentum states are distributed symmetrically about 
the Dirac point (0, 0) with the lowest energy states at (ki, k 2 ) = 

Large gauge transformations in this system correspond to inserting magnetic flux 
quanta (Zi = ^- § dxiAi, I2 = § dx 2 A 2 ) through the two non-contractible curves of the 
torus where the magnetic flux quantum is $0 = 27r/e. The initial choice of l\ and l 2 has 
no effect on the physical spectrum or the allowed momenta: they are all gauge equivalent. 
When l\ or l 2 is adiabatically changed from one value to another, however, we must pass 
through intermediate flux configurations that cannot simply be gauged away. At a generic 
intermediate point in the variation, there is an effect on the boundary conditions for the 
fermions, which now feel a Berry phase due to the magnetic flux as they encircle the 
torus. If the flux through the two holes of the torus is ($i,$2) ? the allowed momenta 

To see the anomaly in action, let us track the fermion spectrum under a large 
gauge transformation in which we simultaneously increase the flux through the two non- 
contractible curves from to 2n. It is necessary to insert non-trivial flux through both 
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cycles; otherwise, there is no zero crossing. As we increase the flux from (0,0) to (tt, 7r), 
the lowest energy fermion states at (kijkz) = (±-^-,±-^-) are shifted so that they sit at 
(ki,k 2 ) = (ffjff 1 )) (0, f^), (f|>0), and (0,0), as shown in Fig. la. If we now continue 
slowly inserting flux, the energies of the states at (|^, |^), (0, |-^), (j^, 0) will remain neg- 
ative, since they are separated by a non-zero gap from the positive energy states in the 
band above. However, the energy of the state at (0,0), which is exactly 0, will change 
sign, so that this state arrives at its final momentum (-^-, -^-) with positive energy, having 
moved up to the conduction band (Fig. 1 b). Because of the single zero crossing, the 
fermion determinant changes sign. 



(a) ^(o.o) (t>) ^(Tt.O) ( a) Y (0,0) ( b) 




Fig. 1. Flux insertion and Kramers degeneracy for non-interacting fermions on the surface of the torus. 

How do we know that a filled state from the lower band must move to the upper band 
during this process? Let |0) be the many-body ground state living on the spatial torus 
and call the operator that inserts (-tt, 7r) flux through the two cycles of the torus $7^. Con- 
sider the two many-body states <Ev )7I -|fi) and T$ n , n \Q). Kramers degeneracy requires that 
(T^ Kj7T Q\^ nt7T Q) = 0. The reason is that the state transported to the tip of the Dirac cone 
after an insertion of (71, n) flux is orthogonal to its time-reversed partner. (For example, T 
transforms the spin-up state at momentum (0,0) to its spin-down partner at (0,0).) We 
then apply to each many-body state to obtain ^^^^^^Q) and ^ niK T^ n!n \Q) ~ |fi). 
The resulting many-body states remain orthogonal. Since the many-body ground state 
with an integral multiple of 2tt flux inserted through each cycle is non-degenerate, the 
resulting states can be orthogonal only if at least one of the states below the Fermi level 
is now empty. In other words, there has been a zero crossing. 
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More generally, the fact that there are always an odd number of zero crossings during 
the above flux insertion process follows from a 3+ld Atiyah-Singer index theorem BS 



Let r parametrize the "time" over which the flux insertion process occurs. Using the 



auxiliary r direction, we may construct the 3+ld Dirac operator = 7 M (<9 M — ieA 



for (A = 0, 3 with A 3 = and Y = a 3 <g> Y for % = 0, 1, 2 and Y = a\ ® \iy.i- The 
determinant of D 3 is obtained from the square root of the determinant of d\ A \ The 
number of zero crossings under an adiabatic flux insertion process is equal to the number of 
zero modes of the four- dimensional Dirac operator The index theorem ensures 



that the number of zero crossings equals A/i ■ A/2. It is important that the change in flux 
number Ali through each cycle is non-zero; otherwise, there would be no zero crossing. 

3.2. The Chern- Simons Term 

When there is an even number of Dirac cones (Nf is even) , there can be no sign change 
of the fermion determinant under any large gauge transformation. If Nf is odd, however, 
we require a bare half-integral level CS term to compensate for the anomaly. 

To see this, consider £7(1) CS theory at level k on the spatial torus S 1 x S 1 , 



Scs = ^~ I e^pA^Ap. (3.2) 
4tt ./(si)3 



The gauge coupling e has been set to unity in this section. We impose periodic boundary 
conditions x\ = x\ + Li, x<i = X2 + L2 on the gauge field and the condition that A^ is pure 
gauge as t — > ±00 implies that our spacetime is effectively S 1 x S 1 x S 1 . (The time direction 
can be thought of as the unit circle if we require Ap i (x,y,t — > —00) = A^(x,y,t — > 00) 
up to gauge transformations.) Consider a field configuration with l\ flux quanta passing 
through the non-contractible curve along x%: 



1 



dx x A x = 2tt/i, (3.3) 



Let us now adiabatically insert I2 flux quanta through the non-contractible curve in X2- 
To do this we must generate an infinitesimal electric field -E^) = < t < 1 (for 

sufficiently large L2). The change in the Chern-Simons term over the course of this flux 
insertion is 

SS C s = 77- / dt f dx 1 dx 2 A 1 E 2 (t) = 2irhl 2 k. (3.4) 







2tt 

k plays the same role as This shows that if k is a half-integer, Scs changes by tt if we 
insert a single flux quantum through the both non-contractible curves along x\ and x 2 - 
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Thus, both the Chern-Simons contribution to the partition function at k — 1/2 and 
the fermion determinant for Nf odd change sign under large gauge transformations with 
l\ = 1% = 1. To avoid the parity anomaly, the overall sign of the gauge field partition 
function must be invariant under these large gauge transformations. This is the origin of 
the condition ( |1.1|) . 

The relationship between the parity anomaly of a single surface and the time-reversal 
protected surface states of the STI can be established without explicit knowledge of the 
surface band structure, as we discuss in Appendix C. This ensures that it remains valid if 
we tune the system away from vanishing chemical potential, or in the presence of disorder 
or strong interactions. Essentially, the statement that each surface, when viewed as an 
isolated system, violates the parity anomaly constraint (|1 . 1|) is equivalent to the statement 
that a bulk system is a STI. 

The analysis above readily extends to more general (continuous) gauge groups. Indeed, 
the parity anomaly was originally discussed in the context of SU(2) gauge theory where 
the analysis is mathematically somewhat simpler 0,^,01 • Let us briefly review this analysis 
as well since the language and notation will be useful later. 

For a non-abelian gauge theory, say SU(2), we may think of spacetime as being 
topologically the 3-sphere by imposing appropriate boundary conditions. (The anomaly 
cares only about topology and not about whether the metric is Lorentzian or Euclidean.) 
In general, we can think of gauge transformations as maps from spacetime into the gauge 
group. Such maps are characterized by their degree or "winding number." The possible 
winding numbers are determined by the group Il3(SU(2)) = Z and so non-trivial large 
gauge transformations are elements of non-zero degree. (The winding numbers of a U(l) 
gauge transformation were denoted by l\ and ^0 

For SU(2), the parity anomaly arises because of the following two facts. First, if the 
SU{2) gauge field is coupled to Nf flavors of 2-component Dirac fermions, then under 
large gauge transformations of degree n, the fermion determinant transforms by (—1) 



iN 



f ■ 



or equivalently, the gauge field effective action shifts by TinNf. Second, if there is a bare 
Chern-Simons term in the action for our gauge field, this term also shifts by 2-nnk under 



a transformation of degree n ||40|| . The combination of these observations again gives us 

Q. 

n 



3.3. Large Gauge Transformations and Domain Wall Fermions 

The discussion above highlights the importance of the topology of spacetime in deriv- 
ing ( |1 . 1| ) . The essential difference between a topological insulator and these purely 2+ld 
systems is that while the low energy fermions are localized to a boundary or defect in the 
fermion mass (assuming a large bulk gap), the gauge field is free to propagate throughout 
the bulk. (See Appendix A for a review of this fermion localization.) Thus, we must 
consider the role played by the bulk geometry. 

Let us gain intuition by first studying the case where the gauge group is SU (2), where 
the geometries involved are simpler. There are two general cases to consider. First, take 
the topology of spacetime to be the 4-ball. This space has a single 3-sphere boundary. If we 
were to forget about the interior of the ball, we could again classify gauge transformations 
by Us(SU(2)) = Z. It is the gauge transformations of non-zero degree that lead to a 
quantization condition on the CS level k. What becomes of them when we fill in the 3- 
sphere to make the 4-ball? Only transformations of degree zero extend continuously into 
the 4-ball bulk. This means that non-trivial (from the perspective of the boundary theory) 
large gauge transformations are not allowed in the theory as they do not continuously 



extend; therefore, (|1.1|) need not be imposed. 

The second situation occurs when the system has topology, say, S 3 x I, where I is 
the unit interval. In this case, the boundary has two components, each living on one end 
of the interval. A smooth extension is now allowed, however, a gauge transformation of 
degree different from zero acts in precisely the same way on each boundary. Because the 
anomalous transformation associated with each S 3 boundary component transforms the 
path integral at most by a phase, the phases associated with each boundary, being equal 
and opposite (due to their opposite relative orientations), precisely cancel one another. 
Because the cancellation is automatic, there is no non-trivial constraint to impose. 

For a general four-dimensional geometry with some number of boundary components, 
the existence of large gauge transformations follows from a basic result of homotopy called 
the extension ,ennn, I A m ap admits a continuous extension if and o„. y if its tota. 
degree (i.e., the sum of its degrees restricted to each boundary component) vanishes. This 
condition ensures that large gauge transformations (if they are allowed with non-trivial 



6 This lemma is familiar from the construction of Wess-Zumino-Witten (WZW) terms. Such 
a WZW term is constructed by extending a field over an auxiliary higher- dimensional space. 
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degree) have no effect on the partition function of any such theory as any accumulated 
phases (e.g., ( |3.4|) from a CS term) must mutually cancel, i.e., there is no anomaly. 

The application of this logic to a theory with U(l) gauge group is now straightforward. 
In the previous section, we considered a 2+ld theory living on x x S^ y y There are 
two simple ways to form a four-dimensional geometry. First, we may fill in the interior of 
one of the circles of the 3-torus so that the entire system only contains a single boundary. 
Without any loss of generality, we may work in A t = gauge and consider the field 
configuration, A x = 2ttIi/L x and A y = 2irl2t/Ly for < t < 1. This configuration is pure 
gauge (i.e., it can be written as A a = e~ l ^ x ' y ^d a e' l ^ x,y ^) at times t = and t = 1. It 
has non-zero winding around and interpolates in time from an I2 = to an I2 = 1 
winding around Sj y If we choose to fill in Sly then A x is no longer pure gauge if l\ 7^ 
as it cannot be removed by a gauge transformation continuously extendable into the bulk. 
This means that only gauge configurations of zero winding in x should be included in the 
physical theory. Similarly, filling in either the or S} y \ is only compatible with AI2 = 0. 
A quick way of drawing the above conclusion is that, in the absence of any sources, the 
Wilson loop around S} x \ or the flux through x must vanish if any of the 
circles is the boundary of some disc. 

Alternatively, we may consider the four-dimensional geometry (S 1 ) 3 x I. This is 
the geometry we implicitly used in the previous section when we discussed flux insertion 
arguments. Large gauge transformations are allowed, but they impose no constraint on 
the physical theory as their effects on each boundary mutually cancel. 

In summary, QED3 must satisfy (|1 . 1|) by having integral ^- + k, if it is to be invariant 
under large gauge transformations. The gauge field coupled to the low energy degrees of 
freedom of a topological insulator, however, is not localized to the boundary; rather, it 
explores the entire bulk. Therefore, we must determine whether or not non-trivial large 
gauge transformations of the boundary theory extend smoothly into the bulk. If they 
do not, then the constraint imposed by their existence is lifted. In situations for which 
they do extend, the number of boundaries is even (or at least the sum of the "degrees" 
of the maps restricted to each boundary vanishes) and the transformation acts identically 
on each boundary, so that (|1 . 1| ) is trivially satisfied. That is to say, there is no non-trivial 
constraint. 
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4. Regularizing the Boundary Theory 

Regularization provides a complementary way to understand the parity anomaly con- 



straint (1.1). In this section, we provide a fairly detailed treatment of the regularization 
of the leading terms in the one- loop effective action for the U(l) gauge field coupled to 
the fermionic degrees of freedom of a topological insulator. Our goal is to discuss how the 
choice of regulator can in general affect the form of this action. In particular, we show 
that the effective action has half-integral ^f - + k in the limit that the bulk gap mo — > oo. 

In general, such a calculation has both IR and UV divergences. To make the calcu- 
lation well defined in the IR, we endow the localized fermions with a mass that is small 
compared to the bulk gap. The mass acts as an IR regulator and allows us to completely 
integrate out all matter fields. However, the mass breaks P and T, and generally induces 
a CS term on the domain wall at level 1/2 (in the limit that the domain wall mass is 
vanishingly small compared with the bulk gap). It is important to consider whether or 
not this level is modified by the regularization procedure. We show explicitly how the UV 
divergences of the boundary or defect theory can be removed without a shift of the CS 
level, and so there is no parity anomaly. In particular, in the limit of a vanishing P and 
T violating mass for the boundary modes, no CS term is generated by the regularization 
procedure - there is no zero magnetic field Hall conductance. This is to be contrasted with 
the situation of a purely 2+ Id theory. 

Because the section is rather technical, let us provide a brief overview. We begin by 
recalling the situation in 2-fTd where the removal of UV divergences in a gauge-invariant 
manner results in the parity anomaly. After this review, we turn to the calculation of the 
domain wall effective action. This proceeds in two steps. We first calculate the propagator 
for the modes localized to the surface. It turns out that while the propagator behaves in the 
IR as one would expect for a localized Dirac fermion, the UV behavior is softer, decaying 
faster at large momentum. Consequently, its contribution to the effective action remains 
finite as the UV cutoff is taken to infinity. (A more careful justification of this claim can 
be found in Appendix E). The bulk is crucial for this effect since the localized modes 
eventually may mix with the bulk continuum if they are excited to energies comparable 
to the bulk gap. Given this behavior, it is then possible to show that a gauge-invariant 
P and T preserving regularization can be chosen. While the discussion concerning the 
regularization of the theory is technical, it has a clear physical interpretation: when P and 
T are locally broken on the surface, our result makes it clear why a level 1/2 CS term is 
obtained for the U(l) gauge field as opposed to an integral level (in the limit of vanishing 
a/mo). 
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4-1- 2+ld Regularization Review 

Let us now briefly review how the parity anomaly can be understood in the context of 
the regularization of QED3 . We will have occasion to make reference to this calculation 
later. Begin with the action, 



S = J d 3 xx(iT(d a -teA a + a^ X , (4.1) 



where x is a single 2-component Dirac spinor of mass a, 7 a = (a 3 , — za 2 , ia 1 ) are 2+ld 
Dirac matrices. Generally, we expect the regularized effective action for A a at energies less 
than a to be a sum of a Maxwell and CS term. 

The effective action is found by calculating the fermion determinant obtained by inte- 
grating out the fermions. It is sufficient to consider the leading quadratic terms in e (and 
consequently in A as well) in the expansion of the determinant. Working in momentum 
space, 



a^A r*,^. i , t i (42) 



S F (A) =e 2 J ^A a (-q)Il<*(q)M<i) 

7o~l3 tr w a 7 - c , n 7 / , x ~ d 
(2ny V (p c 7 c + a) {p + q)dJ d + & 

The kernel U ab (q) or gauge boson self-energy is UV divergent by power counting. 

Thus, we arrive at the delicate (and technical) question of how we should regular- 
ize H ab . The general prescription is to choose a regularization that preserves as many of 
the symmetries present in (fO]) as possible. The usual choice of dimensional regulariza- 
tion does not result in a gauge-invariant effective action; no light fermions are present to 
compensate for the level 1/2 CS term present in the resulting action. The discussion in 
Sec. 3 demonstrates that the CS level must be integral in a purely 2+ld gapped theory 
in order that ( |4.1|) be invariant under large gauge transformations. We are interested in 
maintaining this invariance in our theory so we must choose a different regulator. 

The next (archaic) choice is Pauli-Villars. Here, we first impose an UV cutoff A on 
the momentum integral determining Il ab to find 

U ab (q) = dArf 6 + S -^l e abc iq c + -{q 2 rf h - q a q b ) + 0(1/A), (4.3) 

where c\,c% are non-zero, finite constants. The UV cutoff can be taken to be inversely 
proportional to an underlying lattice spacing. The first term is the UV divergence which 
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manifests itself in a gauge non-invariant mass for A a . Such a term generally appears when 
a strict cutoff is applied to loop integrals in a gauge theory. (Dimensional regular ization 
simply sets c\ =0.) The second and third terms are a level 1/2 CS term and the Maxwell 
term, respectively. 

Next, we introduce a single Pauli-Villars field coupled to A a . The propagator for 
this field has the same form as that for \ except for the replacement a -H- M with a <C 
M. Additionally, the Pauli-Villars field is taken to have bosonic statistics. Thus, its 



contribution to U ab is identical in structure to (|4.3|) except for an overall sign change of 



all terms. Adding these two contributions together gives, 

n „ t = S gn(.)- S gn(M) +ea M-a _ b + 

87T am 

We can now take the cutoff A — > oo. Notice, however, that decoupling the Pauli-Villars 
field leaves behind a non-zero contribution to the CS level. This "spur" is the manifestation 
of the parity anomaly. Parity (and time-reversal) are anomalous in the sense that if a — > 0, 
the starting action is classically P and T invariant. Maintaining gauge invariance in the 
quantum or regularized theory results in the breaking of P and T if we are to use the same 
regulator for both the massless and massive theories. 

Within this regularization scheme, it is not possible to remove by some clever choice 
of Pauli-Villars fields the gauge non-invariant UV divergence without a non-zero shift of 
the CS term. This conclusion is in complete agreement with the topological argument 
of the previous section and we will, therefore, adopt Pauli-Villars as our regularization 
prescription. That is, the requirement of invariance of the low energy action under large 
gauge transformations is satisfied within the Pauli-Villars regularization scheme, but not 
within dimensional regularization. 

The lesson is that there is a physical difference between dimensional and Pauli-Villars 
regularization in 2+ld. If invariance of the theory under large gauge transformations is to 
be maintained, then Pauli-Villars regularization must be used and parity is broken. From 
the discussion of the previous section, we anticipate that there is no physically observable 
difference between Pauli-Villars and dimensional regularization when the theory lives on 
the boundary of a contractible higher-dimensional space because there do not exist large 
gauge transformations that extend continuously into the bulk in any such theory. In the 
next section, we show explicitly how this equivalence comes about. 
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4-2. Propagator 

We now turn to the calculation of the regularized effective action for the localized 
modes described at low energies by ( [4 . 1| ) . Lest we run into a contradiction with the 
conclusion of Sec. 3, the action ( [4.1] ) must only be a low energy approximation to the 
physics. We will find that the 3+ld bulk, which is implicit in any low energy domain wall 
action, will play an essential role. 

In order to derive the action, we first determine the full propagator for the domain wall 



modes by generalizing the nice work of Chandrasekharan |[41|| . Perhaps surprisingly, the 
form of the propagator is partially determined by the massive bulk modes. Our derivation 
is contrasted with an anomaly inflow argument in Appendix B. We have also repeated in 
Appendix D the calculations in this section for the analogous system in the simpler 1+ld 
context. 

As we reviewed in Sec. 2, we describe the 3+ld topological insulator by the action, 



S 



d 4 x ^[i^^d^ - ieA^) + m(x) + i^a^ip, (4.5) 



where ijj is a 4-component spinor that describes the massive bulk fermion with spatially 
dependent mass m(x) coupled to the U{1) gauge field which may represent the elec- 
tromagnetic field. The spatially-varying mass allows us to study the interface between a 
topologically trivial and topologically non-trivial insulator, while the constant T-breaking 
mass a regularizes the infrared divergences on the domain wall. (We have not included 
kinetic terms for the gauge field, since their specific form does not affect our results.) Note 
that tp = -i/^7 . 

We take m{x^) to depend only on a single coordinate x^ and to have the profile, 
m(x3) = m o tanh(x3/£) with mo > 0. If a 7^ 0, we may safely integrate out %j) without an 
IR divergence. We assume that < a <C mo- (In the opposite limit, there is no domain 
wall and the entire system is in a single massive P and T breaking phase). 

We are interested in the one-loop effective action for A^ obtained after integrating 
out ijj. This action takes the generic form, 

S = ~y J d^xd A x'A lJL {x)A v {x')tx^D{x,x') 1 v D{x',x)^ (4.6) 

where D(x, x') is the fermion propagator and the overall minus sign comes from the fermion 



loop. We are particularly interested both in divergent terms which are implicit in (|4.6| ) that 



require regularization and in possible CS terms localized to the 2+ld domain wall. These 
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CS terms may arise either directly from modes localized to the domain wall or through an 
integration by parts of a bulk #-term. 

Thus, our first task is to calculate D(x,x'). We choose the following representation 
for the Dirac matrices, 

"_ ;/ M J"' <M .,5 _ ,.,0.4.,2.,3 /<> I 



V = . 7' = "i (^o _;, J ,7 s = JtVtV = (J J . (4.7) 
The eigenspinors of the Dirac equation at e = 0, 

7°(*7 M ^ + m(x 3 ) +z 7 5 a)v A = A^ A . (4.8) 
formally define the propagator, 

Q(*y) = iV^ W f V) , (4.9) 

A A 

where the eigenfunctions ifj x are normalized with respect to the usual inner product, 

d 4 x(ij x )\x)ij x '(x) = 5(X - A')- (4.10) 



Translation invariance in the temporal and spatial directions parallel to the domain 
wall allows us to express the propagator in a mixed Fourier space representation, 

D(k, X3 , X ' 3 )=ij: ^f^\ (4.1D 

A 

where k collectively refers to momentum parallel to the wall and ifrfc are the Fourier coeffi- 
cients in the expansion of i/j x . Expanding in Fourier modes in this mixed basis, the Dirac 
equation becomes, 

l°(l a k a + z 7 3 d 3 + m(x 3 ) + za 7 5 )^ = A^, (4.12) 

where a = 0, 1, 2. Two types of eigenstates are expected: modes localized near the domain 
wall and modes allowed to propagate away from the domain wall. We refer to the former 
as bound modes and the latter set as scattering modes. 

Consider first the scattering modes. The eigenvalues for the scattering states may be 
obtained by considering (|4.12| ) in the limit £3 — > ±00. In this limit, translation invariance 
along the X3-direction is effectively restored and so we introduce an asymptotic momentum 
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k 3 parametrizing the eigenvalues of Q4.12| ). There are two pairs of eigenspinors ^^ 3 -' ± 'W 
with eigenvalues, X(k 3 )± = X± = ko ± \/k\ + k\ + k\ + m\ + a 2 = ko ± u k - 

The coupled first order equation ( |4. 1 2| ) may be rewritten as the second order equation, 



(dl + kt + ml(l + ^)sech 2 (^))^ = 0. 



(4.13) 



Above, we have substituted the mass profile m(x 3 ) = ttiq tanh(x 3 /£). This equation says 
that each component of tp^ satisfies a modified Poschl- Teller equation, with parameter 
mo/£. The associated eigenspinors and eigenvalues to this equation are known for general 
mo/£, however, they can be given simple closed form expressions when I = I/vtlq. For this 
reason, we specialize to this point in parameter space for the remainder of the paper. 

There is a well known connection between the Poschl- Teller equation and supersym- 
metric quantum mechanics (see e.g., |4^]). The upshot is that the separation of eigenspinors 
into bound and scattering states can be made precise. Eigenspinor components of scat- 
tering states are paired; while those of the bound states are zero modes of the associated 
supersymmetric charge operator. The number of such zero modes depends upon the ratio 
£/mo which we set to unity, thereby implying a single set of bound state modes. See 
Appendix D for further details about this connection in the technically simpler case of 
1+ld. 

An orthonormal basis for the scattering state solutions to ( [4.13| ) is provided by 



/ k 3 - im(x 3 ) \ 



V 



ki + ik 2 
±u k - a 




(kt - ik 2 )(k 3 - im(x 3 )) \ 



(4.14) 



-kl - ml 



^±,(2), , __L 

(±u k - a)(k 3 - im(x 3 )) I 



where 



A/l 1} = y/2(2ir)*u k (u k Tcr), Af^ = Af^ J k% + m 2 . 



(4.15) 



It remains to find the bound states. A bound state takes the form ip^° 



axi bx2 ) sech(moX3). These states do not carry "momentum" k 3 ; they are 
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localized to the domain wall because of the prefactor sech(moa;3). A normalized basis is 
provided by 



ki - ik2\ 

^° ,±iX3) = ^ 4(2,)C>T g ) ( . 

±t k - a/ 



sech(moX3), (4-16) 



with eigenvalues ko ± tk = ko ± \/k\ + k\ + a 2 . 

Having found the eigenspinors and eigenvalues, we can now construct the propagator 
using (|4.9|) . It is a sum of two terms arising from the bound and scattering modes, 



D(k, x 3 , 4) = L> bound (fc, x 3 , 4) + D scat (k, x 3 , 4). (4.17) 
We compactly write the contribution to the propagator from the scattering states as 

Z^,*3,4) = i [ M ( XS,X> 3 ) +i*l 5 eik3(x >- X3) (41g) 

J 2-k k* — fcg — ttiq — a 2 
for fx = 0, 3 and where k% = Kq — k\ — k\. The unconventional mass matrix, 

M(x 3 ,x 3 ) = — (I+17 ) g - (l-*7 )+2( 1+z 7 ) (fc2 + m 2) M^,^)- ( 4 - 19 ) 

and 

M^3, 4)fc 3 = ^(^3)^(4) + ik 3 (m(x 3 ) - m(x 3 )j - mg. (4.20) 

Notice that M(a;3,4) approaches the standard form —mo if the mass m(x 3 ) is taken to 
be a constant. The bound state contributes the following term to the propagator, 



.bound/ t\ im n ^-^ W \ W ,^ a ka + H 5 ° 

I {x, x ) = - -^-(1 + «7 )sech(m a;3)sech(moa;3) — % — 



for a = 0, 1,2. 



4-3. Corrections to the Effective Action 

We now use the above propagator to study corrections to the gauge field effective 



action. Instead of plugging the full propagator into (|4.6| ), it is useful to first study the 
structure of the propagator itself more closely. 

Before doing so, two technical comments are in order. First, we have been working in 
Lorentzian signature in the previous sections, however, we have found it most convenient 
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to perform the necessary (intermediate step) integrals below by first Wick rotating to Eu- 
clidean signature by substituting ko — > iko so that k 2 — > — \k 2 \. Our resulting expressions, 
however, are written with the original Lorentzian metric. Second, we use a renormalization 
procedure in which we integrate over all k^, but impose an UV cutoff A on the remaining 
three momentum integrals k a . This unconventional choice should not affect the low energy 
properties of the regularized theory. 

inbound is the Fourier transform of the usual 2+ld Dirac propagator. Indeed, de- 
spite the slightly different matrix structure, substituting ]j hound into (fO|) gives precisely 
the same contribution to the effective action as a truly 2+ld Dirac fermion reviewed in 
Q!.2|) . The hyperbolic prefactors merely localize the contribution to the domain wall. The 
localization becomes exact in the tuq — > oo limit where we use the relation 

lim — ^-sech 2 (77102:3) = S(xs). (4.22) 

mo-K» 2 

Therefore, we must understand how the UV divergences are regularized. In particular, we 
wish to determine whether or not there is a shift to the CS level if we regularize using a 
Pauli-Villars scheme. 

A clue comes from a closer inspection of the full propagator. Not only do the bound 
modes result in propagation that is localized to the domain wall; surprisingly, the scattering 
modes contribute a term that is localized as well! 

The unconventional terms in the mass matrix M(xs,x' 3 ) proportional to ii(xs,x' 3 ) 
are responsible for this localization. Let us momentarily focus upon these terms in the 
scattering mode contribution Q4.18|) to the propagator proportional to /i(x3, x' 3 ). We denote 
these terms by _D sc . at After a contour integration over k^, they take the form 



e ik 3 (x' 3 -x 3 } 



(4.23) 

where k^ above is evaluated at k% = zsgn(x 3 — xs)^/— k\ + m 2 , + a 2 . 

Remarkably, the scattering modes contribute a term that is equal and opposite to the 
contribution to the propagator arising from the bound modes. These two terms cancel 
one another in the full propagator. The remaining localized term in the propagator is the 
second term appearing in (|4.23| ). Thus, we define the localized propagator, 

fl-(*,**4) = -!^) (1+i7 3) x k : + ,r\ s^"-', (4.24) 

4 (kl - a 2 ) y/-kl + ml + a 2 
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where, as above, k 3 = isgn(x' 3 — x 3 ) a/— k 2 + wig + a 2 . This propagator describes the 
excitations that are restricted to living along the domain wall at energies low compared 
to the bulk gap. This is the main technical result of our paper and it is the form of this 
localized propagator that is the key allowing us to understand the divergence structure of 
the 2+ld boundary theory. 

Notice that the UV behavior of D loc is softer than £) bound : it decays faster as k — > oo. 
This behavior implies that the purely local mode contribution to the effective action ( [4.6| ) 
is finite as the UV cutoff is taken to infinity. (By purely local mode contribution to the 
effective action, we mean (|4.6|) with D loc substituted in place for the general propagator 
D.) We shall begin by calculating these finite terms. Afterwards, we comment upon the 
corrections to the action arising from interactions between modes localized to the domain 
wall and those modes that are free to propagate throughout the bulk. A detailed analysis 
of these latter two terms is relegated to Appendix E. 

To isolate the finite contributions to the action from the localized modes, we find it 
convenient to use the following approximate expression for the local propagator, 

D a P pToAk,x 3l x 3 ) = ™(l + i7 )— sech(m a;3)sech(moa;3) 



2 V ' 2 v v (^- ( t 2 )V-^ + ^o + ^ 2 ' 

(4.25) 

This expression becomes a better approximation to the exact result ( [4.23| ) as the limit 
m — > oo is approached. We expect corrections to this approximation to be suppressed in 
the large bulk mass limit. (In Appendix E, we evaluate the contribution of the localized 
mode without making the above approximation and find that ( f4.25| ) gives the correct 
qualitative structure for the divergences and the quantitatively correct value for any finite 
terms). 

Using the approximate local propagator ( |4.25| ), we may now compute its contribution 
to the gauge field effective action. In particular, the kernel or one-loop self-energy is 

U ab (q, x 3 x> 3 ) = ( - l -^rj cd + + ^(qV 6 - Q a q b ) + —))8(x 3 )S(x' 3 ), 

V 6 87T \a\ a z mo J 

(4.26) 

where C2 is a finite non-zero constant and the delta functions arise from the hyperbolic 
prefactors using ( f4.22| ). We stress that the above result is finite due to the faster decay 
of (|4.25| ) at large momentum, so that a 2+1-dimensional UV regulator is not required. 
However, we have exchanged a gauge non-invariant term proportional to the cutoff for 
precisely the same term, now proportional to the bulk gap mo. 



22 



Therefore, we must understand how to properly decouple the bulk by taking uiq — > oo. 
It is essential that the divergence is proportional to the bulk gap rather than an UV cutoff. 
This difference allows us to regularize the theory by Pauli-Villars fields without shifting 
the CS level. (In Appendix E, we show in detail that the linear divergence appearing 
as mo — > oo is the only possible divergence that we may associate with presence of the 
localized states. 0(1) corrections to the coefficient of this term from the value displayed in 
( f4.26| ) may arise from interactions with the bulk modes, however, the precise value of the 



coefficient does not affect the arguments below. Additionally, there are the usual (bulk) 
divergences associated to QED4.) 

We merely need to introduce two Pauli-Villars fields whose propagators have roughly 
the same form as the localized fermion (|4 .25 ). The differences lie in the choice of statistics 



Ci and masses for the Pauli-Villars fields. We endow the first Pauli-Villars field with 
bosonic statistics, C\ = 1, and replace m B- M and a B- a' . The second field is taken to 
have fermionic statistics, C<i = —1, and we replace mo B- M — mo and a B- a'. We assume 
the hierarchies uiq <C M and a a' < mo- Adding the contributions from the physical 
boundary fermion and the two Pauli-Villars fields, ( 4.26|) becomes 



rr t M 1 / I i.^i^i.^i^ | sgn(q) + Cisgn(q / ) + C 2 sgn(a') , cbd 
Uab^x^Xz) = -(- two + Ci\M\ + C 2 \M - m | H iq b e , 

(4.27) 

where we have suppressed writing the contributions to the Maxwell term and it is to be 
understood that the above correction is localized at X3 = x' 3 = 0. Clearly our choice 
of statistics and masses removes the term in ( |4.26| ) proportional to mo, but retains the 
half-integral level of the CS term induced by the localized fermion. As promised, we have 
regularized the boundary theory with Pauli-Villars fields in a P and T invariant manner, 
i.e., we have regularized without shifting the CS level and so no CS term is generated by 
the regularization procedure when a — > 0. We were successful in doing so because of the 
softer UV properties of the localized fermion which resulted in a divergence proportional 
to the bulk gap rather than the UV cutoff. This enabled the theory to be regularized with 
an even instead of an odd number of Pauli-Villars fields. 

The form taken by the Pauli-Villars fields' propagators betrays their 3+ld origin. 



Indeed 3+ld QED requires at least three Pauli-Villars fields for its regularization [g3[. We 
have merely chosen two of the Pauli-Villars fields in the 3+ld bulk to have soliton masses 
similar in form to the physical fermion so that they can regularize the lower-dimensional 
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theory. The remaining 3+ld Pauli-Villars fields can be given masses that are positive 
everywhere. 

The finite gauge non-invariant photon mass term (which we removed in ( 4.271) ) may 
be surprising. However, its appearance is similar to what occurs in the Pauli-Villars 



regularization of QED4 where a gauge non-invariant mass term is also found [43[] . In that 
case, the mass squared is a sum of two terms: one proportional to the square of the cutoff 
and a finite term proportional to the square of the bulk fermion mass. In both cases, the 
appearance of gauge non-invariant terms proportional to a power of the fermion mass is 
a result of the fact that imposing a momentum cutoff A (in our prescription for the 3- 
momenta k 2 , or for all 4- momenta k 2 in the standard QED4 case) breaks gauge invariance 
explicitly. The correct choice of Pauli-Villars regulator fields is determined by the criterion 
that they must fully restore the gauge symmetry broken by this choice of cutoff. 

Now that we have explained how the excitations described by the local propagator 
are regularized, we should ask: Can the massive bulk modes contribute non-trivial terms 
to the gauge field effective action? By massive bulk modes, we mean the terms in _D scat 
that are not localized to the domain wall by any hyperbolic prefactors. Specifically, the 
full fermion propagator is given by D = D loc + _D free , with D loc given by ( f4.24j ), and 



(4.28) 

where as above, k% = i sign(#3 — x^)-\/—k 2 + m\ + a 2 . The corrections that we have thus 
far ignored arise from either single insertions of D loc and _D free or two insertions of _D free 
into ( f4.6| ). These correct both the boundary and bulk Lagrangians. For example, the bulk 
Maxwell term is radiatively modified by the massive bulk modes described by D free , while 
the cross-term of D loc with _D free modifies both the 2 + Id Maxwell term and the effective 
Chern-Simons term arising at the domain wall. 

As with the local contribution discussed above, there are two aspects to this question. 
First, we must determine if there can be a direct correction to either a bulk #-term or (what 
is the same) a boundary CS term. Second, we must determine what effect (if any) these 
contributions have on the choice of regulator fields that we must include in the theory. 

It is straightforward to show that the only correction to the boundary CS term arises 
from the "crossterm" between D loc and D free . No bulk #-term is generated. Adding this 
crossterm correction to our result ( [4.27| ), we find the total CS level, 



2\a\ n \nrno 



), (4.29) 
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where r\ = ±1 is defined by the orientation of the domain wall, m(xz) = r\m§ tanh(mo£3). 
(Our calculation above specialized to the case r\ = 1. Similar manipulations show that 
Qf.29| ) obtains for rj = —1.) Precisely the same non-quantized correction to the CS level 
also occurs in the analogous 1+ld situation - see (D.23). (In 1+ld, the CS level directly 
determines the induced charge on the soliton, however, in 3+ld the soliton only carries a 
charge when a background magnetic field is applied perpendicular to the surface.) Extrap- 
olating the 1+ld intuition (axion electrodynamics is suggestive as well |y]|), we expect 
that for general (not necessarily constant) a, 

k = f* (4.30) 

where A(j) is the change in phase <fi of the complexified bulk fermion mass during any 
interpolation. The fact that ( |4. 3Q|) need not be rational is not in conflict with gauge 
invariance; for a single domain wall, there is no quantization condition on the CS level (as 
shown in Sec. 3). For a system with two domain walls, the phase of the fermion mass must 
wind by an integer multiple of 2-k as it crosses both boundaries, provided that the fermion 
mass asymptotes to the same value everywhere outside the system. 

Are there divergences appearing as the cutoff A — > oo that must be regularized? We 
study these possible divergences in Appendix E and show that they do not occur. There 
are two types of possible divergences. The first kind arise from the usual massive 3+ld 
propagator and are familiar from QED4. Such divergences are always present and can 
be regularized (with at least three Pauli-Villars fields ||43|1 ) without changing the effective 
CS level. The second type of divergence is special to the introduction of a domain wall 
mass profile and is localized along the wall. Such divergences are at worst logarithmic in 
A; however, a precise cancellation occurs between the bulk and localized modes, described 
above, such that only finite terms remain as A — > oo. We stress that this cancellation 
is important: a divergent result would have required adding an odd number of regulator 
fields with domain-wall mass profiles, thereby rendering the CS level k an integer (up to 
the corrections of order a /mo noted above), forcing the constraint ( |1 . 1| ) to be obeyed in 
the limit a /mo — > 0. 
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5. Conclusion 

In this work, we have examined the relationship between the parity anomaly and the 
gapless Dirac fermion arising at the surface of a 3+ld topological insulator. Naively, a 
coupling of the localized surface mode to a fluctuating gauge field would result in the parity 
anomaly on each boundary surface, i.e., the requirement that + k be integral on each 
boundary surface where Nf is the number of fermions localized at the boundary, and k is 
the sum of any bare or induced CS level. We have shown that this does not occur and that 
^y- + k is half-integral in the limit that the bulk gap mo — > oo. Hence, an odd number 
of gapless Dirac fermions on any such surface can be coupled to fluctuating (bulk) gauge 
fields and still maintain parity and time-reversal symmetries. 

We have come to this conclusion using two complementary perspectives: topological 
quantization conditions and regularization. It is invariance of a 2+ld effective theory 
under large gauge transformations that results in the constraint that + k be integral. 
When the system of interest lives on the boundary of a higher-dimensional space, such 
large gauge transformations either do not exist (as they cannot be extended continuously 
into the bulk) or they are innocuous - their effect is cancelled between all components of 
the boundary. Thus the topological quantization conditions that gauge invariance imposes 
on 2+ld theories do not apply in 3+ld. Similarly, Pauli-Villars regularization of QED3 
preserves the invariance of the theory under large gauge transformation at the cost of 
breaking parity and time-reversal invariance. We have shown explicitly how for a single 
species of Dirac fermion on the 2+ld boundary of a 3+ld bulk, the presence of the bulk 
softens the UV properties of the boundary fermions, such that the theory admits a parity 
and time-reversal invariant Pauli-Villars regularization. 

It is instructive to contrast our result with the quantum Hall effect. Here, invariance 
under local gauge transformations (charge conservation) requires that both a bulk CS term 
and boundary chiral excitations be present in a low energy description of the system. In 
contrast, only a global anomaly can be present for the surface modes of a 3+ld topological 
insulator. However, the fact that the gauge field is free to explore the bulk essentially 
eliminates the large gauge transformations responsible for the possible global anomaly; 
there is no mutual cancellation of anomalous transformations as occurs in the quantum Hall 
effect. Thus unlike the chiral edge modes of a quantum Hall system, the gapless boundary 
modes of a 3+ld topological insulator are not required to preserve gauge invariance. 

The fact that there are no topological quantization conditions for these boundary the- 
ories has important physical consequences. A purely 2+ld system exhibiting a fractional 
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Hall conductivity of a xy = -\ with integer p can be gauge invariant only if the system 
exhibits a ground state degeneracy equal to p when placed on a spatial torus. Other- 
wise, the theory fails to be invariant under large gauge transformations generated by flux 
insertions through the two non-trivial cycles of the torus. If we fill in the center of the 
torus, however, the non-trivial large gauge transformations that cause this problem no 
longer exist (in the sense that they do not continuously extend into the bulk), and the 
requirement of a degenerate ground state disappears. Hence, as we would expect for a 
non-interacting system, the ground state of a topological band insulator is unique in spite 
of the fact that p = 2 on its surface (in the limit of vanishing a /mo with finite a). (For a 
fractional topological insulator in 3+ld, there is a ground state degeneracy, but this arises 
from its bulk topological order [|23j, and not from the surface). If instead we thicken the 
torus, gauge invariance requires a ground state degeneracy that is half of what one might 
naively expect based on the Hall conductivity, since the large gauge transformations must 
behave identically on both surfaces. For the topological band insulator this again implies 
a unique ground state. 

In fact, the perturbative calculation reveals that when both bulk and boundary 
fermionic mode contributions are considered, the coefficient of the Chern-Simons term 
(physically, the Hall conductivity at the surface) is not quantized to be a rational fraction 
when — is non-vanishing. Rather, it is of the form 

mo ° ' 

fc^-Itan- (-£-), (5.1) 
2\<j\ n \rjmQ/ 

where rj = ±1 determines the orientation of the domain wall via the soliton mass profile, 
m(x) = ??mo tanh(moX3). Here, mo is the bulk band gap, and a is the time-reversal 
and parity breaking mass of the boundary fermions. For a general soliton mass profile, we 
expect the CS level to equal l/2it multiplied by the total change in phase of the complexified 
fermion mass. (In Appendix D, we show that the CS level for an analogous 1+ld system 
is also given by ( |5.1| ).) We may think of the half-integral contribution coming from the 
fermionic modes localized to the boundary, while the non-quantized contribution arising 
from the interaction between the bulk fermions and the modes localized to the boundary. 

In a real system where the ratio a /mo is not asymptotically vanishing, the above 
deviation from half-integrality may be observable. A magnetic-susceptibility measurement 
would in-principle measure the difference in the Hall conductivities. Combined with a 
Kerr or Faraday rotation measurement (which measures the sum), one could then extract 
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the Hall conductivities of each boundary surface [33]. However, we caution that strict use 
of the formula (|5.1|) requires a non-vanishing time-reversal breaking perturbation to be 
present at asymptotically large distances from any topological insulator boundary. 

We stress that in order for our conclusions to be valid, the gauge field must be free to 
explore the bulk of the system. It is possible to imagine a situation in which strong corre- 
lations among the modes localized to the boundary of the system lead to fractionalization. 
If fractionalization only occurs at the boundary, it can be described by an "emergent" 
gauge field that only has support on the lower-dimensional boundary and so any resulting 
constraints imposed by gauge invariance truly are of a lower-dimensional origin. An alter- 
native scenario for localizing a truly 2+ld gauge field uses a bulk Higgs field charged under 
some gauge group whose symmetry-breaking profile only allows a subgroup (possibly, an 
empty one) of light "photons" in the bulk, but a domain wall defect where the full gauge 
symmetry is restored. We hope to discuss these scenarios further in future work. 
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Appendix A. Domain Wall Fermions 

In this appendix, we review the domain wall fermions P3],H that arise in the model, 



S = J cfx^^O^ -ieA^)ip + m(x)ipipJ, (A.l) 
where the mass m(x) is real and satisfies, 

lim m(xz) = ±mo, rn > (A.2) 

X3— >-±oo 



7 We thank S. Kachru for discussions on this latter possibility. 
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with m(x 3 ) passing through zero exactly once, at x 3 = 0. We can think of these gapless 
domain wall fermions as zero modes bound to a defect in the order parameter where 
symmetry is restored. In this example, the order parameter is the real field m(x 3 ) and 
the symmetry that is (classically) restored at the defect or location where m{x 3 ) vanishes 
is the chiral symmetry, ijj — > exp(ia n /5)ijj . (This is only a true symmetry of the quantum 
theory at zero gauge coupling due to the chiral anomaly ]T3,IB|,0].) 



To see how these zero modes arise, it is useful to write the Dirac equation as, 

(iD\\+iD ± +m(x 3 )\ip = Ey ip, (A.3) 



where 

D\\ =^ a {d a -ieA a ), a = 0,1, 2 
D ± =~f 3 (d 3 -ieA 3 ). 



(A.4) 



A zero mode satisfies the equation, 

(iD ± +m(x 3 )^il) = E~f°il), (A.5) 

with E = 0. We are interested in finding a solution about the free theory so we may set 
e = in the equations of motion. 

It is convenient at this point to choose a particular representation for the Dirac ma- 
trices. The results, however, are independent of any particular choice. We take 

l2 <r)^=(-°< <m> 

where a 1 for i = 1, 2, 3, are the usual Pauli-sigma matrices. Note that the above choice is 
a rotation from that used in Sec. 4. 

Since tr(7 3 ) = and (7 s ) 2 = — 1-4x4) 7 3 possesses two pairs of eigenvalues equal to 
±z with eigenspinors defined by the equation 

7 3 ^± = ±itp±. (A.7) 



A solution to (|A.7|) contains half the number of degrees of freedom of a 3+ld Dirac fermion. 



The eigenspinors take the form 

±ia 6 x 
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^±(x)=<P±L\ ), (A. 



where x( x ) is an arbitrary two-component spinor that depends only on the three coordi- 
nates parallel to the domain wall and 

P ^3 

<f>±(x 3 ) = C±</4 0) = C±exp (± / m(x 3 )). (A.9) 

C± is a normalization constant and x^ can be chosen to coincide with the location of the 
domain wall. In the following analysis, it is convenient to normalize ip± by choosing C± 2 = 
2 J (0^ ) 2 dx3. While both ip± solve the zero mode equation, only one is normalizable; only 
one has a finite, non-zero C±. The asymptotics ( A. 2 ) chosen above for m(xs) singles out 



ifj_ as the normalizable zero mode. Had the opposite asymptotics been chosen, would 
be the normalizable zero mode. 

Now consider the action that describes these fermionic localized modes at energies 
much less than the bulk band gap, ttiq. It is found by substituting i(j± into the 3+ld 
action ( |A.1|) . Given the profile (|A.2|) for m(xs), we plug ?/>_ into the action to find 



--C 2 / dz (6^) 2 



S(ip-) = J d x[ip-i'f^(d f j i — ieA /Jr )ip_ + m(xs)ip-4'- 

rf 3 x((x t -^(a 3 x) t )7 7 a (9a-^eA a )(x ia 3 xY 
= J d*x(x H a (d a ~ ieA a ) x y 

(A.10) 

The above action describes a massless 2-component Dirac fermion coupled to a U(l) gauge 
field, namely QED3. Again, we have suppressed the tree-level kinetic term for restricted 
to the domain wall. It is sufficient to say that the tree-level gauge boson propagator 
restricted to the 2+ld surface is softer in the IR, decaying as l/\p\ as opposed to 1/p 2 at 
small momentum, because of an integration over the direction normal to the domain wall. 
The resulting 2+ld Dirac matrices are 

7 a = ( a 3 , -ia 2 , ia 1 ) , a = 1, 2, 3. (A.ll) 

Note that the coefficient of the minimal coupling term in the action between the fermion 
number current in the 23-direction and A3 vanishes identically since 

^±7 3 ^± = 0. This 

ensures that no zero modes leak off the defect via a coupling to A3. 
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Appendix B. Anomaly Inflow 

In this appendix, we first review anomaly inflow for the case of a string defect in 3+ld. 
We then contrast this analysis to that of a domain wall defect in 3+ld which is relevant 
to this paper. 

The original model studied by Witten [|45| and Callan and Harvey |14[ 



is 



S = J d^x ift(i^{d„ - ieA^) + m(x)e l ^ip. (B.l) 

where ip is coupled to a string defect defined by the complex field m(x) exp(i(ft). m(x) 
vanishes at the core of the string running along the z-axis and (ft winds by 2n in going 
around the string. This mass profile ensures that tere exist chiral fermionic zero modes 
living on the string which exhibit a gauge anomaly through their coupling to A^. 



However, integrating out the massive bulk modes provides a Wess-Zumino term ||46 
whose gauge variation compensates for the lack of gauge invariance of the chiral zero 
modes alone. This can be understood through the following calculation. Consider a region 
of spacetime away from the string core where m(x) is non-zero and let us integrate out ift. 
The correction to the action by the massive bulk fermions can be inferred by integrating 
the one-loop correction to the current expectation value, 

g 

(Ju) = pW 9 ^^' ( B - 2 ) 
Because (ft winds by 2n when encircling the string, it formally obeys the equation, 

(d x d y - d y d x )(ft = 2n5(x)5(y). (B.3) 

Thus, (|B.2|) implies the bulk fermion addition to the current conservation equation, 

= ^F tz S(x)5(y). (B.4) 



The current is conserved away from the string and the right hand side of (|B.4j ) is non- 



zero along the string so as to cancel the contribution from the chiral zero modes. This is 
summarized by a correction to the effective action, 

e 2 f 

Swz = / d A xe^ vpa d^(ftA v F pa , (B.5) 
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whose variation under a gauge transformation cancels with the anomalous variation of the 
chiral zero mode action. The cancellation between boundary and bulk anomalies is called 
anomaly inflow as the direction of current flow is towards the boundary. 

We stress that the contribution of the massive bulk states is added to the low energy 
effective action for the string defect. We shall not find this prescription to be strictly valid 
for a domain wall defect. 

In order to describe a domain wall, we merely change the soliton profile. m(x) now 
vanishes along the domain wall and <fi jumps from to n in moving through the wall. 
Naively, precisely the same calculation of integrating out the fermions in a region where 
m(x) is non- vanishing gives the contribution to the effective action ( P-5| ). 

However, this reasoning is in fact incorrect. The essential difference is that <p is con- 
stant everywhere except in a small neighborhood of the domain wall where m(x) vanishes. 
This is to be contrasted with the string case where <p wound by 2tt around the string. 
Because </> is constant in the region of space where the calculation obtaining (p.5| ) is valid, 
this correction vanishes everywhere that the calculation is well-defined. 

We can make the calculation well-defined in the vicinity of the domain wall if we 
smooth out the field <p. Instead of jumping discontinuously at the location of the domain 
wall, we allow it to smoothly interpolate between zero and 7r. In effect, this smooth 
interpolation imparts a non-zero P and T breaking mass that we denoted by a in the main 
text. As we have understood through more direct means, this merely breaks P and T on 
the defect and results in a level 1/2 CS term on the domain wall. 

The discontinuous limit is not strictly available within the above scheme. This limit 
is equivalent to taking a — > 0. Indeed, integrating out the fermions is not well defined in 
this limit when the mass of the localized fermions vanishes and so the a — > limit need 
not commute with this integration. Instead, if we are interested in vanishing bound state 
mass, we first take a — > and then integrate out the fermions to obtain the 1PI effective 
action. This action does not contain a CS term. 

Appendix C. The Parity Anomaly and Strong Topological Insulators 

In this appendix, we briefly recall the relationship between the presence of time- 
reversal protected gapless surface states and the parity anomaly in strong topological 
insulators. 
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This relationship is most apparent from the definition of a STI proposed by [28,23|. 
They consider the fate of a STI on the thickened spatial torus S 1 x S 1 x i", which has 
two disconnected toroidal surface boundaries, and two non-contractible curves through 
which we may insert magnetic flux. There are four flux choices for which the system 
is time-reversal invariant ($i,$2) = (0, 0), (tv, 0), (0, 7r), and (7r,7r). If the many-body 
ground state is Kramers degenerate (i.e., orthogonal to its time-reversed conjugate, which 
is necessarily a state of the same energy) in an odd number of the 4 flux sectors, the surface 
spectrum has 2+1 dimensional gapless surface states (in the thermodynamic limit) [B3|. 
These gapless states cannot be eliminated without breaking the Kramers degeneracy, and 
therefore breaking time-reversal symmetry; hence the system is a STI. This definition is 
equivalent to the band-structure based definition of [fHHI in the non-interacting case, but 
has the advantage that it does not require an explicit knowledge of the band structure, 
and thus is equally applicable to interacting systems. 

This criterion ensures that, if time-reversal symmetry is preserved, there must be 
a large gauge transformation in which at each surface an odd number of fermions cross 
from below the Fermi surface to above it. Let us begin in a flux sector where the ground 
state is non-degenerate. Next, we adiabatically insert flux to arrive in a flux sector where 
the many-body ground state is Kramers degenerate. Applying T (which maps the many- 
body ground state to its orthogonal Kramers partner), and inserting the same flux again 
must return the system to its original many-body ground state (since this is equivalent to 
inserting no flux at all) . If we simply apply the flux insertion twice, without performing a 
time-reversal transformation in between, we therefore obtain a fermionic configuration that 
is orthogonal to the original. As the original many-body ground state was non-degenerate, 
after this large gauge transformation the system must be in an excited state. If the only 
zero-energy state encountered during the flux insertion is in the Kramers doublet, then 
only one fermionic mode crosses the Fermi surface. (More generally, extra band crossings 
not protected by Kramers theorem can occur in pairs at momenta (p,—p), and an odd 
number of fermions can be transferred.) Hence, there exists a large gauge transformation 
in which an odd number of fermions cross the Fermi surface. 

Appendix D. 1+ld Effective Action Calculation 

In this section, we repeat the calculation of Sec. 4 in 1+ld. This is technically simpler 
than the higher-dimensional calculation and so is easier to follow. Our results are consistent 
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with those in Sec. 3 and also allows us to make direct contact with the work of Goldstone 



and Wilczek 26 



We begin with the 1+ld action, 

S = J d 2 x tp^i^id^ -ieAp) + m(x) + ij 5 (rjip, (D.l) 

where tp = ip^l°i 0" is a constant, and lim x ^± 00 = ±mo for uiq > with a single zero 
crossing at x = 0. Any such soliton configuration in the mass is stable in 1+ld as opposed 
to higher dimensions. If a ^ 0, we may safely integrate out if) without any IR divergences. 
This produces an effective action for at energies less than a. The leading term is linear 
in and we expect it to take the form, 

S eS = -e J d 2 x ^tr^Dfox)), (D.2) 

where D(x,x') is the propagator for ifi evaluated at x = x' (note that D(x,x') is a 2 x 2 
matrix and so the trace does not automatically vanish). Without resorting to an argument 
26|,|14j presented in Appendix B, it is necessary to construct the real space propagator 
because there is no translation invariance in the direction perpendicular to the domain 
wall. 

Thus, we must calculate D(x,x'). First, we choose the following representation for 
the Dirac matrices, 

/ = a \ 7 1 = «t 3 j 7 5 = 7 V=<t 2 . (D.3) 

The propagator of the Dirac fermion is determined by the eigenvalues and the eigenfunc- 
tions of the equation, 

l°(il^d„ + m(x) + z 7 5 o)^ A = A^ A , (D.4) 

which can be written as 

idt-a d x +m(x)\ f ip x \ _ ^ fl>i\ ( D5 ) 



— d x + m(x) idt + a J \ip2 J \' l / J 2 
The propagator is formally given by the expression 

ip x (x)ip x (x') 



A 

A 

where the ip x are normalized eigenfunctions. 
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D(x,x')=ij: ^ ^ (D.6) 



Time-translation invariance allows us to Fourier expand t^f (t, x) = f ifti(w, t) exp(-iwt). 
Thus, we must solve the following eigenvalue problem, 

-d x + m(x) u + a j \^ A (w,i) / \ "02 ( u > x ) J 

Spatial translation invariance is restored as x — > ±00. We can use this fact to solve for 
the eigenvalues (of any states that are not localized to the domain wall) by asymptotically 
Fourier decomposing, linr^-i-oo ^(u, x) = f, ^(w, k) exp(— ikx). (p.7| ) becomes 

u-a -ik±m \ (ip$(u,k)\ = x (ip£(u),k)\ (J) . 

ik±m u + a J \ip$ (uj , k) J \ip^(u,k) J ' ^ '' 

Thus, there are two sets of eigenfunctions, , ± (a;,x), with eigenvalues, A± = ui ± 
a//c 2 + m\ + a 2 = uj ± ojk- In addition to these scattering states, there are modes bound 
to the domain wall as we shall describe. 

To find the eigenfunctions of the scattering states, we must solve the equation, 

u - a d x + m(x) \ ( i>i(u, x) \ _ . ( ipf(u,x) 



o x + m[x) ui + a J \ip 2 (uj,x) J \tp 2 {00, x) J 

It is possible to do this by finding the eigenfunctions for the constant mass Dirac equation 
and then replacing the constant mass by the spatially varying one. However, the following 
observation is instructive. The coupled first-order equation (|D.9| ) may be written as two 
decoupled second order equations, 

d 2 — d x m(x) — m(x) 2 + m\ + k 2 )i/jf =0 

' (D.10) 



(d 2 + d x m(x) - m{x) 2 + ml + k 2 ^jipf =0. 



To proceed further, we must choose a particular form of the interpolating mass. We take 
m(x) = mo tanh(x/£). Substituting this into ( P-lOj ) we find 



(^ + fc2 + (1 - 1/TOo ,)_^L)^ = 

2 (D.ll) 

(d 2 x + k 2 + (l + 1/moi) ™° - W =0- 
V cosh lx/r.)s 



cosh (x/£) 

These are generalized Poschl- Teller equations. 
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There is a well known relation between the 
ric quantum mechanics as we now describe [[12 
Define the (supersymmetric charge) operator, 



Poschl- Teller equations and supersymmet- 
. For simplicity, let us work at I = I/tjiq. 



Q = ds — m tanh(m a;3), Q' = —83 — mo tanh(m £3). 



(D.12) 



We may suggestively rewrite (|D.11| ) so that the equation satisfied by each eigenspinor 
component takes one of two forms: 



QiQWfh = ("2 - ^)(^ A± )i, QQW) 2 = («J - a 2 )(V fc A± ) 2 , (D.13) 



where the subscript, 1,2, refers to the spinor component. Recall that the eigenvalue in the 
above equation {uf, — a 2 ) = k 2 + m\. 

Scattering states are defined to be those eigenf unctions for which (uf. — a 2 ) ^ 0. 
There is a pairing between scattering state spinor components; given a solution to the 



identical eigenvalue is given by Q4>i- A bound state is annihilated by either Q or Q< and 
there is no corresponding pairing of spinor components. Consequently, bound states can 
only become scattering states in pairs. The number of bound states is given by the index of 
Q; this number is the difference in dimensions of the kernel (or null space) of the operators 
Q and Q f . For our problem, only states annihilated by Q* are normalizable and so we 
need only consider this operator when finding the bound state wave function. The number 
of bound states or zero modes is preserved under small changes of £. For < m§l < 1, 
there exists a single zero mode bound state. 

Following this brief digression, let us now directly solve ( p.ll| ) at the point 1/i = uiq 
for which the two equations simplify to 



The equation for t(; 1 ± is easily solved by const. exp(— ikx). It may be substituted back into 
the first order differential equations ( p.9|) in order to determine ip^ '• Doing so, we find 
the un-normalized eigenspinors for the scattering states, 



first equation in (|D.13| ) with non-zero eigenvalue, a solution to the second equation with 




(D.14) 



±u k - a 
ik + rug tanh(mox) 



) 



—iwt—ikx 



(D.15) 
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Note that k only has the strict interpretation of momentum asymptotically far away from 
the domain wall when m(x) becomes a constant; otherwise, it may be viewed as a repa- 
rameterization of the eigenvalues. It is interesting that a never explicitly enters the second 
order differential equations (|D.11|) . 

It remains to find the bound state. The two components of the Dirac spinor are no 
longer paired. Instead, 1^2 is the only non-zero component of the spinor and we need only 
solve the first order differential equation ([D . 9| ) with m(x) = mo tanh(mox). (There is no 
corresponding normalizable solution where ^2 = and ifii is non-zero.) The un-normalized 
solution is 

^ (^)=(sech('L ,)) e " i " i ' ( D ' 16 > 

with eigenvalue Xq = u + a. It is a bit of a misnomer to refer to ifj x ° as a zero mode as its 
energy is bounded from below by a. It is better to call it a bound state. Of course, this 
bound state is free to mix with the bulk continuum when its energy is comparable to the 
bulk gap. 

Equipped with the eigenspinors (p.!5| ) and ( p.!6|) , we must now normalize them with 
respect to the inner product ( |4.10|) where we replace the number of integration dimensions 
by two and n = 1, 2 depending upon whether the state is localized or extended. Thus, the 
orthonormal eigenspinors are (same notation as un-normalized ones above) 



An V sech(moa;) 



U^uJkiukTo) V^ + m tanh(m ^ > 6 



Let us now construct the propagator using (p.6| ) . This is a 2 x 2 matrix. The sum over 
the eigenvalues becomes integrals over u, k. The zero mode only contributes a non-zero 
21-entry (given our choice of 7- matrices), 

D hound (x,x') = ^(l + ^7 1 )sech(mox)sech(mox , ) [ < ^- 1 \ [ + ^ e^' ~ l \ (D.18) 
4 J 2tv uj 1 — a 2 

Note that this is simply the propagator of a massive 0+ld particle. 

The contribution to the propagator from the scattering states is more complicated, 
but straightforwardly found as before, 

)S cat, A / dujdk 7 M ^M + M ( X ' X ') + Au(t'-t)+ik(x'-x) 



D scat (x,x') = / , , 2 ; ^ ?2 2 2 eM } X ~ X ' ( D - 19 ) 
' (2-7r) 2 uj z — k 1 — ttiq — a 1 
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M(x, x') = -hl- z 7 1 )m(x / ) - i(l + iT 1 )™^) + ^(1 + H 1 ) "?"*^ /^ x ') ( D - 2 °) 

^ ^ A K "i 

and 



fj,(x,x') = m(x)m(x') — ik(m(x) — m(V)) — mQ. (D.21) 

The similarity between the two-dimensional and four-dimensional cases is evident. 

Return to the expression for the leading contribution to the effective action for A^, 

Seff = -e J d 2 x AfJx^Dfax)). (D.22) 

In 1+ld, we only need the propagator evaluated at the same starting and ending point. 
Therefore, let us simply let t' = t and x' = x in our expressions for the bound and scattering 
state propagators. Using our expression for the propagator at coincident points, we find 

S M = - Itan-'(^)) fdt eA , (D.23) 

where r\ = ±1 allows for a general domain wall orientation, m(x) = r/rrio tanh(mox). We 
also used the relation, 

, m sech 2 (m x) B , N _ 

lim — = S( x ). (D.24) 

mo— too 2 

( p.23| ) is simply the 0+ld CS term for the Aq field. In the limit of interest, a /mo — > 0, the 
level is half-integral. Notice that the total level depends upon the asymptotic ratio a /mo 
precisely in the way predicted by Goldstone and Wilczek ]2H[] . While our calculation was 
performed for a particular soliton profile, we expect the level to be l/2n times the change 
in phase of the bulk fermion mass. 



Appendix E. Bulk and Boundary Corrections 

Here we will give the details of the calculation of the 1-loop correction to the gauge 
field effective action. This correction takes the form, 



SS =J d 4 xd 4 x'A^x)W v {x,x')A v {x') 



e 2 



J d 4 xd 4 x'A^(x)A l/ (x')tr(^y^D(x 1 x , ) r y v D(x' , x) j . 



(e.i; 



In Sec. 4, we studied in detail the contribution of D loc to ( |bC.l| ). Because of the softer 
UV behavior of this propagator, we found this term to be finite in the UV cutoff. Thus 
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the level 1/2 CS term arising from the fermion propagator is, unlike in a 2 + Id system, 
unmodified by any regularization procedure. 

It is important to verify that any divergences in the A — > oo limit arising from other 



contributions to (|E.1|) do not force us to add regulator fields that will change the level of 
the Chern-Simons coefficient modulo 1. To verify this, we write the 3+ld propagator as 
the sum of two terms, 

D(k, a*, 4) = D loc (k, s 3 , 4) + D iiee (k, s 3 , 4), (E.2) 

where 

(E.3) 

and 

n free/L / \ _ l a kg-^K-\ [m(x 3 ) + m(x' 3 ) + z 7 3 (m(s 3 ) -m(s 3 ))] + za 75 _ fc3 | a ./_ 3 . 3 | 

±> (ft,S 3 ,S 3 J - — F== 2,2 6 

2 V + mg + a 2 

(E-4) 

for fc 3 = \f—k 2 + rriQ + a 2 , k = ik 3 sign(x 3 — S3). Here we have integrated over all k 3 
in ( }4.18| ), such that the contribution from the poles at k 3 = ±irriQ exactly cancels the 
propagator of the bound states, as explained in Sect. 4. This gives n Ml/ = nj^+nj^+n^ 
with 

11$ = Tr ( 7M D free 7 „D free ) = Tr ( 7/J D loc 7 „D loc ) (E.5) 

Tl% = Tr ( 7M D free 7l ,D loc + 7M D loc 7 ,D free ) . (E.6) 

In the main text we discussed only IT],,; here we will evaluate the remaining contributions. 

To understand the structure of the possible divergences, it is helpful to Fourier trans- 
form in S3. Recall that the assumed translation invariance in the directions tangent to the 
domain wall allow us to immediately work in a mixed Fourier basis where k a for a = 0, 1, 2 
is the momentum tangent to the domain wall. Focusing on the third spatial direction, we 
introduce the momenta s, t: 



dx 3 dx' 3 Yl^(x 3 , s 3 , (f a )A M (s 3 , q a )A„(x' 3 , -q a ) 
= / dx 3 dx 3 / ds dt e lsXa e ltx ^A^(s,k a )A v (t,k a )Ii llv (x 3 ,x' 3 , k a ) 
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Intuitively, this form is convenient because fixing x 3 , x' 3 at the location of the domain 
wall requires us to include modes of arbitrarily high energy, obscuring the structure of 
the divergences. Thus, it is preferable to first integrate over x 3 ,x 3 and then examine the 
structure of the divergences for small s,£, q a . 

As a sanity check, it is instructive to see what this prescription gives for the case of a 
translationally-invariant mass. In this case, the only x 3 dependence of is 

g — (^3+^3)1^3— 053 1 
C3/C3 

where k = (k , fei, k 2 , k 3 ), k! = (k - q , h - qi, k 2 - q 2 , k 3 ) and 



/_ e -(fc3+fc 3 )|x3-X 3 | 
d 3 k a -j—j- f ((k ■ k' - m 2 )r]^ - k^k' u - k'^k u ) , 



fc 3 = \J-k 2 + k 2 + k 2 + m 2 + a 2 , k' 3 = ^-{k - q ) 2 + (h - Ql ) 2 + (k 2 - q 2 ) 2 + m 2 + a 2 . 

(E.8) 

In this case, performing the integral over the center-of-mass co-ordinate x 3 + x 3 enforces 
the condition s = —t. The remaining integral over x 3 — £3 gives: 

J d(x 3 -x 3 )U^(x 3 ,x 3 ,p a ) = J d s k a 2 ^ k + >y + s 2 ~j~^r ' k ' ~ m ' 2 ) 1 l^ - KK ~ k '^ k ») • 

(E.9) 

As one might expect, this is exactly the result we would obtain by first Fourier transforming 
the propagator in all four spacetime momenta, and then integrating over k 3 , with q 3 = s. 

When the mass breaks translation invariance in x 3 , s and t are not conserved, and 
integrating over x 3 , x 3 will not force s = —t. However, it remains true that U /Jil ,(x 3 , x' 3 , q a ) 
is an admixture of modes at many different energies, and that we must integrate over x 3 
and x 3 in order to obtain the correct structure of divergences. 

Defining z = x 3 — x 3 , Z = x 3 + x' 3 , the various contributions to the propagator are: 
r e -(k 3 +k' 3 )\z\ 

n ^("?a, x 3 , x 3 ) = / d 3 k a j—jp [(k p (k p + q p ) - ml - a 2 )r]^ - fc /i (fe I/ + q v ) - K(k^ + g M ) 

1 . „ \ 2Tf 



+ VTo - T\ ^ + i {k3 + k , 3) sign(z) ^9a5,3 ) - S, 3 S v3 - {h + k , y 
r e -(k 3 +k' 3 )\z\ 

^l(Qa,x 3 ,x 3 ) = J d 3 k a -j-p T cross [5^ a 5 ub {(k c (k c + q c ) - a 2 )r] a)b - k b (k a + q a ) 

-k a (k b + q b ) + iae ajbjC q c } - 5 p3 5 u3 (k c (k c + q c ) - a 2 ) + a] 

f e -(k 3 +k' 3 )\z\ 

11^(90,^3,^3) = J d 3 k a -j—^, T loc [S fia 5 ub {(k c (k c + q c ) - a 2 )r] a , b - k b (k a + q a ) 

-k a (k b + q b ) + iae a ,b,c<lc}] 

(E.10) 
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where we have defined k' 3 = y— (fco + Qo) 2 + (k± + q±) 2 + (fa + Q2) 2 + m o + c 2 - Here a is 
a matrix of linear order in q a that is non-zero only in the fourth row and column, and 

mf. . 2 ,2 / \ o 2 1 + cosh mo Z cosh uiq z 

T Q = — (tanh m x 3 + tanh m x 3 ) + m = 2m - — — 

2 (cosh moZ + cosh moZ) A 

m (k 3 + k' 3 )sign(z) , m (k 3 + k' 3 ) smhm \z\ 

±1 = (tanhmox 3 — tanh mox 3 ) — 



2 " cosh mo Z + cosh mo z 

2 / 1 1 / \ 2tyiq cosh z /-rp 

T% =m (l — tanh moX3 tanh moXo) = l-^- 11 ; 

cosh mo Z + cosh mo z 

2~cross = , 



, 2 

1 



For k 2 large, 11^ and 11^ scale like while 11^ scales like \k\~ 2 - hence naively, 

after integrating over the remaining three loop momenta, all three might be divergent. 
However, care must be taken with this naive power-counting, as is apparent from the form 
of n^ D4 prior to integrating over z, Z: power-counting suggests that the leading-order 
divergence should be cubic, while it is in fact quadratic in A. Likewise, we will find that 
the leading-order divergence from the terms that arise due to the spatial variation of m(x 3 ) 
is logarithmic. 

To exhibit these divergences explicitly, we next Fourier transform in x 3 , as in Eq. 
( |F.7| ), and perform the integrals over z and Z. For example, we wish to integrate 

i J dZ dze-^We^+W^e^-W'To. (E.12) 

To do this for general s and t, we observe that the zeros in the denominator of To occur 
at Z = ±z + i(2n + 1)tt. Summing over n, we obtain 



- I dZ e ~ 2 k 3 \z\ i(s+t)/2Z i(s-t)/2z T 

2 J 



1 _ e -2k 3 \z\ e i( S -t)/2z 



(E.13) 

[Res(T ; x = y + in) + Res(T ; x = y - in)] . 

We then integrate these expressions with respect to z. After performing both integrations, 
the mass dimension of the result has been reduced by 2 (rather than by 1, as it is when 
integrating over Z leads to S(s — £)). This indicates that the corresponding contribution 
to n M „ has the mass dimension of the correction for a 2 + ID theory. 
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in 



We obtain: 

dz dZe- 2ka \ z \e i{s+t)/2Z e i{s - t)/2z T = 27rsign( 



m J 



ks(s + t)(kl + Us 2 +t 2 )) 



(E.14) 

The exponential e^^" 1 "*-*/ 7710 ensures that we may expend the result for small not only 
do these parametrize the photon momentum in the X3 direction, which we take to be small 
relative to mo, but equally the result is exponentially suppressed in while after Wick 
rotating, k% > ■sj Wq + a 2 . 

This expansion can be safely performed for all of the integrals. To zeroeth order in 
s, t, we obtain: 



dzdZe- 2k ^Tn 



2|mo| 



i J dzdZe- 2k ^T x 



\m \ 



[ dzdZe~ 2k ^\T 2 = f-^f + 2\m Q \d k3 r(-^ 
J \ k 3 I mo I 



dzdZe 



- I dzdZe 



-2k 3 \z\rr2 



ml + 2/c 3 |m | - 2/c3|mo|9fc 3 r(- 



\m \ 



-2k 3 \z\rp2 



Tj=2(2k£ + m£ 



2/c| 


m 


+ 1k\ml - 8fc| 


mo 


3 + k 2 rriQ + 4/c 3 mo 


5 - 4m 3 


H(k 3 -2\ 


>™ol) 


2 (h- 


m |) 2 



fa 



+2|m |a fc3 r(-2 + ^)) 



I w 1 



(E.15) 

where T is the digamma function. 

Taylor expanding these expressions for large k a , one can see that nj]},(0, s = 0, t = 0) 
is non-divergent, as claimed in Sect. 4. Further, the divergent terms associated with the 
domain wall in IT^(0, s = 0, t = 0) and IT^(0, s = 0, £ = 0) cancel. 

We can also evaluate the extra finite contributions to the gauge non-invariant terms 
due to the contributions we neglected in Sect. 4. Inserting the expressions in ( E.15Q 
into the expression for the Fourier-transformed propagator, and Wick rotating to send 
k 2 — > —k 2 , we obtain: 

rV(0, 0, 0) = IV(0, 0, 0f ED *- [ rf3fcc 2mo(mg + 3ay) ^ = _2jmo[ 5 

J 6K 6 (k~ + o i y 3 

(E.16) 

where IV (0, 0, 0) Q ™ is the gauge non-invariant contribution from the fermion loop in 
QED4. 
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Thus the gauge non-invariant terms that arise due to the domain wall do indeed have 
the form given in Eq. (|4.26| ), and regularization can be carried out with an even number 



of Pauli-Villars fields with spatially varying masses. 

We can also use this method to calculate the coefficient of the Chern-Simons term in 
Eq. ( E.10|) . As described in Sect. 4, U 1 ^ contributes a half-integral Chern-Simons term, 
whose sign is determined by the sign of the domain-wall mass a. However, there is also a 
non-quantized contribution from Ilj^,. Thus the total effective Chern-Simons term is 

n S? = A ( 7rsi § n ( (J ) - 2 sin" 1 a | (ie abcP c ) (E.17) 
or a Chern-Simons coefficient of 

k = -sign(a) - - tan -1 (E.18) 
2 7T \mo\ 

where we have used sin -1 , = = tan -1 t- 2 -,. This is reminiscent of the corrections 

^/ml+a 2 \ m o\ 

to the induced charge of solitons in the simpler 1+ld model ( p.23|) : in 3+ld if we break 



time-reversal (for example by applying a magnetic field to the system) the domain wall 
carries an induced charge, whose value is exactly quantized only in the limit a /mo — > 0. 
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